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1 Introduction 



The analogy between the 2-d XY model and its 3-d counterpart in the sense that both models 
possess a phase transition at finite critical coupling [3^, driven by vortices and anti- vortices, may one 
lead to suspect that the behavior of the vortices and anti- vortices in the 3-d XY model is analogous to 
the 2-d XY model. 

The generally accepted picture of the phase transition in the 2-d XY model is that for couplings 
/3~/3c, the vortices and anti- vortices are present with finite density, and behave effectively as free 
particles. As the coupling approaches the critical coupling, the vortices and anti-vortices start to form 
bound pairs, and may disappear from the system by annihilation, thereby decreasing the vortex density, 
or vorticity. For couplings (3^ (3c, the system is described by spin waves, while vortices and anti- vortices 
are only present as tightly bound pairs at (very) low density [1]. 

If we impose periodic boundary conditions, we will find that in the 3-d case the vortices and anti- 
vortices will have to form closed loops. The importance of vortex loops on the phase transition of the 
3-d XY model has been noted by [2]. The intuitive picture of the phase transition in the 3-d XY model 
is that for couplings [3 > [3^, there will be few small loops present in the system. As the coupling [3 is 
lowered, several possibilities arise. Either the number of loops will increase while the length of the loops 
remains (roughly) constant, or the length of the loops will increase, while this time the number of loops 
reamins constant. Also an increase in both is possible. 

We will study this picture of the phase transition by measuring properties of these vortex loops. A 
preliminary account of this work is given in [3] . 

This paper is divided into sections as follows. In section 2, we will describe the 3-d XY model, 
and introduce the vortex loops. After discussing some of their properties, we will define the relevant 
vortex loop correlation functions and distributions. Section 3 will be devoted to the presentation and 
discussion of the results. In section 4, we will recapitulate these results and give a refined description 
of the phase transition. 



2 The 3-d XY Model 

The action for the 3-d XY model is given by 
S' = /3^ cos((?,+^ - (?,). (1) 

The quantity 

kfivix) = — |[6'j;_|_^ — 9x] + [9x + (j. + v — (^x + li] + [(^x + v — (^x + (j. + v] + [(^x — f^x + v]^, 1^ < t^, (2) 

where [. . .] denotes the restriction to the interval < — tt, tt], yields the vortex number on the corresponding 
plaquette Pfj,v{x). To study the properties of these vortices, we perform a duality transformation on k 
by defining [4] 

i;i(^ ~ /*) = l^^livp kvp{x), (3) 

where x = x + ^fi + + ^p. 

The links on the dual lattice have to satisfy the following relation: 

d'^j,{x) = V£, (4) 

d'^ /(«) = /(«) - /(« - f^)- 

This means that every point x on the dual lattice is visited by an even number of non-zero links j. On 
a finite lattice, this implies that we find closed loops, or clusters of closed loops. 
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Consider now the following spin configuration on an elementary cube: 



■'X — "x + 2 — ^' 

- — - — -ivr 

9 - - — - - — ivr 

. . — _3_ 
-^x + 1+3 ~ 4 ' 

a ... - 3 
-'x + 1+2+3 — 4 



(5) 



See also figure 1. This configuration has a vortex A; = 1 on the right face, and an anti- vortex k = —I 
on the top face. The e-tensor conspires in such a way as to have the same sign for both non-zero fc's, 
leading to, dropping the tilde from now on: 

jl{x) = l, j2{x)=0, j3{x) = -l, 

- i) = j2(x - 2) = js(x - 3) = 0. (6) 

So we see that a single loop may consist of links j = 1, as well as links j = —1. Moreover, a link j = I 
may be related to a vortex or an anti- vortex. If we now assign a direction to the links j, we see that 
the conservation law (4) implies that the number of links pointing to x is equal to the number of links 
pointing away from x. Although for non-intersecting planar loops, we can distinguish between loops 
going around in a clockwise fashion from loops going around in an anti-clockwise fashion, this becomes 
impossible when these loops become non-planar, or intersect a number of times. Therefore, we shall 
make no distinction whatsoever between different loops. 

In order to study the behavior of these loops, we define the following two correlation functions [3,5]: 

Cs{r)= J2 ^-eu^yeu/ E 1, (7^) 

\x-y\=r \x-y\=r 

Q«= E ^^eiAei,/ El- (^E^-e'.)'- 

x,y,i,j x,i 
\x-y\=r \x-y\=r 

Cs(r) yields the probability that two points x and y, a distance r away from one another, are on the 
same loop li. C'^(r) measures correlations between possibly different loops. C's also has a field theoretic 
connection [5]. If we describe these links j by a complex scalar field ip, so jfi(x) = ip(x)ip'^ (x + fi), then 
we find 

Cs(r) = (^+(x)^(x)^+(y)^(y) 



C 



y\ 



where the subscript C stands for 'connected part'. From the fact that we can only have an even number 
of links at each point x we can deduce that only four point and six point interactions are allowed. 

Now, when we have to deal with percolating loops, i.e. loops wrapping around the lattice, we will 
find that C's(r) will tend to a constant for large separations r. This constant is of course related to the 
density of these percolating loops. 

However, for Cj(r), this constant is cancelled by the second term in (7b). So, in general we expect 
the correlation functions (7a) and (7b) to behave like 

^/ X , exp (-ar) , , 

^('■) = d + c. (7c) 

If Cs is dominated by the first terms of the expansion of < Lp'^ {x)Lp{x)Lp'^ {y)Lp{y) >, i.e. eye-like 
diagrams, see figure 2a, we expect from the dimensionality of the field Lp the parameter h to be equal to 
2. However, if also (very) high order terms have a significant contribution to < Lp'^ {x)Lp{x)Lp'^ {y)Lp{y) > 



3 



in this expansion, see figure 2b, tlien it bears more resemblence to a bound state, and tlierefore we 
expect h to be equal to 1. 

In tlie spin wave pliase, we expect tlie constant c in (7c) to be equal to zero. Then, if we do 
not find any difference in behavior between C's and Cj, so in other words, if we find the same values 
for the parameters a and h, then both correlation functions measure the same type of correlations - 
self-correlations. This then implies that the loops do not interact with one another. 

In the vortex phase, we do expect different behavior when the vortex loops percolate, as then C's 
will tend to a constant for large separations, while this constant is absent for Cj. 

As there is nothing to prevent loops from intersecting, we expect loops to do so. However, there is no 
clear way to distinguish between a self-intersecting loop and two loops touching one another. There are 
several ways to deal with this problem. We feel that in such a case, we should average over all different 
ways of dividing up such a loop. The reason is the following. If we start following a vortex loop we will 
arrive at a point where two or three links are pointing outwards. We then have to decide which link to 
follow. This choice will determine whether we follow the entire loop, or split the loop in two (or more) 
smaller loops. This is most easily visualized by considering the shape 8. Arriving at the crossing we 
have to decide in which direction to continue. One choice will lead to following the whole 8, the other 
will lead to splitting the shape up into two small rings. 

Having defined loops, we can also determine the distribution of loop lengths p(l). Again, when 
loops percolate, we expect loops of any length in the system, hence the distribution p(l) will tend to a 
(non-zero) constant for large lengths /. 

3 Results 

We used a 10 hit Metropolis algorithm with site dependent step size, combined with 7 = 2 multigrid, 
as described in detail in [7]. We made runs consisting of 10 batches of 100 measurements each, after 
discarding 500 measurements for thermalization. The measurements were taken after each time the 
spins on the original lattice were updated, i.e., after each half W-cycle. Both C's and were obtained 
from averaging over the 10 batches. 

We made runs for /3 = 0.00 up to /3 = 1.00 in steps of A/3 = 0.05, as well as runs with A/3 = 0.01 for 
/3 between 0.40 and 0.50. Two additional runs were made at /3 = 0.525 and at /3 = 0.575. The lattice 
sizes we used were 6^, 8^, 12^, and 16^. 

To define a vortex loop, we choose the following strategy. We start at some point xq, which happens 
to have a link j z/z 0. This point xq will obviously be on a vortex loop. We then scan in a prefixed 
order for the next point of the vortex loop, marking the connecting link as read. This search is repeated 
until no more unread links are to be found. We have then returned to the starting point xq- We then 
search for a next point xi, and repeat the procedure. If we visit some point x' , to which four (or six) 
non-zero links are connected, there is the possibility that this self-intersecting loop will be split into two 
(or more) different loops, depending on the geometry of the original, self-intersecting loop. Since the 
system is isotropic, there is no preferred geometry, and this procedure will in effect average over the 
different possibilities of splitting up loops (including no splits at all). 

Apart from the correlation functions C's and and the length distribution p(l), we also measured 
the average number of loops the average loop length A, and the fractal dimension df. The latter is 
defined as the ratio of the total number of links over the total number of points visited by the vortex 
loops [5]. 

3.1 Cs and Q 

In figure 3, we show the behavior of Cs on a 16^ lattice. For /3 in the range /3 = 0.40 - 0.50. As can be 
seen from this figure, there is a difference in behavior between [3 < [3^ and [3 > [3^ (Pc ~ 0.4542 [8-10]). 
Notice that Cs approaches a constant for (3'^j3c, but that for [3 > /3c, this constant is absent. When 
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making a x^"fit of Cs(r) to (7c), we find the results as given in tables 1 to 4. It is interesting to note 
that for [3 < j3c the parameters a, h and d are relatively independent of [3, while only c does depend on 
[3. As a matter of fact, c is monotonically decreasing for increasing [3. This means, since c is related to 
the density of percolating loops, that also this density is decreasing for increasing [3, and vanishes when 
c vanishes. 

Figures 4a and 4b show Cs and for a 16^ lattice at /3 = 0.30 and j3 = 0.55, respectively. It is 
observed that falls off exponentially on either side of the phase transition, and that for [3 > [3^, both 
Cs and C(j fall off with the same exponent. We did not fit to (7c) since there were not enough data 
points to make a reasonable fit. 

In figure 5, we have shown the exponent a versus [3 in the range [3 = 0.35 io (3 = 0.60, for all four 
lattice sizes. It can be seen from this figure that a starts to rise approximately linearly for (3 > [3^, up 
to /3 0.60. Figure 6a shows the constant c versus [3 for [3 = 0.00 io (3 = 0.60, while 6b shows c in the 
range [3 = 0.35 io (3 = 0.50. We see that the results for both a and c for the lattice sizes N = 12 and 
# = 16 are within one another's error bars. 

The upshot of these results is that there is no clear sign for interaction between loops on either side of 
the phase transition. In the spin wave phase, we see that Cs(r) behaves identical to Cj(r), while in the 
vortex phase Cj(r) falls off with a large exponent. An estimate from figure 4a gives Cj(r) exp(— 15 r). 

Focussing on Cs, we see that the mass a vanishes in the vortex phase. From the fact that h 1.5, 
we can deduce that the field Lp does form a bound state, albeit a weakly bound one. Finding a non-zero 
value for c in the vortex phase, and a zero value in the spin wave phase means that the system undergoes 
symmetry breaking, and the non-zero value of c indicates that a condensate has been formed. 

The fact that the results for N = 12 are equal to the ones from N = IQ implies that we do not need 
to study larger lattices if we are only interested in the behavior of vortex loops. 

3.2 #, A, and d^. 

In tables 5 and 6, we give the results for A, and df. The errors were obtained by blocking. This 
information is also displayed in figures 7a-d. Several interesting observations can be made. The first 
one is that the density of loops, SI = zf^/N^, falls off exponentially for (3^(3^. This means that loops 
become exponentially suppressed in this coupling region. This is hardly surprising as this region is the 
spin wave phase, where the dominant degrees of freedom are spin waves. Note that here SI does not 
show much finite size effect, in contrast to the vortex phase, [3 <^ [3^, where indeed there are finite size 
effects. The second observation is the similarity between A and df. It appears that, on average, the 
larger the loops, the more they (self)intersect. Again note that the results for N = 12 and N = IQ are at 
most a few a away from one another in the coupling region 0.30 < [3 < 0.60, supporting the observation 
made earlier that these lattice sizes are large enough to obtain reliable results for vortex loops. 

3.3 p{l) and I 

Figures 8a and 8b depict the behavior of p{l) versus / for 16^ lattices at /3 = 0.30 and j3 = 0.60, 
respectively. Again we see different behavior on either side of the phase transition. The presence of a 
constant for /3 < /3c is clear, while for [3 [3^, there is an exponential suppression of larger loops. Notice 
that there still are a small number of loops with lengths a few times the linear extension of the lattice. 
In tables 7 and 8, we present the results for fitting p{l) to (7c). For 16^ we have plotted a and h from 
this fit in figure 9. Note that both a and h behave (roughly) linearly for 0.45 < (3 < 0.60. An estimate 
would give a 2.36/3 — 1.14, and h Ri 6.8 — 9.8/3. The fact that the parameter a becomes non-zero for 
[3 > j3c is, of course, related to the exponential suppression of vortex loops found in the previous section. 

Turning our attention now to the vortex phase, we see that the parameter h is constant, albeit with 
large errors. An estimate of & in this region would give h Ri 2.2(3). For a, we would estimate a Ri 0.01(7). 
These values are consistent with the results from [11]. The results for the spin wave phase do not agree, 
mainly due to the absence of an exponential factor in [11]. 
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As we have made 10 batches per run, we can average over the length of the largest loop per batch to 
get an estimate for the length of the largest loop for a given /3- value. This average length, I, is presented 
in table 9. We have also made a log-plot of these data in figure 10. It should be kept in mind when 
reading tables 7 and 8 that / runs essentially from to I. Notice that I scales differently on either side 
of the phase transition. 

4 Conclusions and discussions 

The picture of the phase transition is now quite clear. As /3 \ /3c, we see that the loop density SI 
increases, while the average length A remains constant, or increases slightly. As /3 > /3c, but quite close 
to /3c, the loop density remains constant, while the average length increases sharply. If we make an 
estimate in this region of how many loops of length A can fit in, say, a 16^ lattice, we find that at 
/3 = 0.48 there are on average 125 loops in the system. The average length at that /3-value is 10.6. We 
can accommodate a loop of such a length in a 3^ block. So we see that about the maximum number of 
loops are present in the system. As [3 is further decreased into the vortex regime, we see that the loop 
density remains roughly constant, while the average length increases roughly linearly with decreasing 
[3. So around the phase transition, the length of the loops is becoming so large that different loops 
start touching one another, and will fuse to form one, large loop. This loop will necessarily have many 
self-intersections, which can be witnessed from the behavior of df in this region, and from the similarity 
between A and df. At the phase transition, finally, these long loops start percolating through the lattice. 
This is substantiated by the fact that for /3 < /3c , C's(r) c ^ Q for large r. 

All this also indicates that the different loops interact rather feebly with one another, if at all. If 
they would have interacted more strongly with one another, we would expect only a handful of larger 
loops at /3- values slightly above /3c, instead of roughly 100 short loops, filling the whole lattice. Another 
clue is the fact that Cs and fall off with the same exponent for [3 > [3^. Figure 4b shows both at 
/3 = 0.55. For this /3-value, there are on average 80 loops in the system. Had the loops interacted with 
one another, we would have seen a difference in behavior between the two correlation functions. This is, 
of course, not very surprising. Since a single loop may consist of links j = 1, as well as of links j = —1, 
we expect no significant interaction between links of opposite sign, nor of like sign. Any interactions 
between links would lead to the collapse of the loops they constitute. 

The onset of percolation also explains the absence of an exponential in Cs and p for [3 < [3^. As 
percolating loops destroy any long range order, these exponentials cannot be accustomed for, and are 
absent, or small, due to finite size effects. 

The description in terms of the field Lp is the following. For couplings above /3c, the system exhibits 
a finite mass a. At /3c, the system undergoes a symmetry breaking phase transition, at which the mass 
vanishes. Furthermore, a condensate is formed, which is witnessed by the constant c. In this condensate, 
the Lp fields form weakly bound states. 

So the picture of the phase transition in the 3-d XY model is the following. At large couplings [3, 
the system is dominated by spin waves, and allows only a few, small vortex loops. As the coupling [3 is 
lowered, the number of loops increases, but the average length remains constant. As [3 approaches /3c, 
the number of loops becomes so large that they start to fuse, and form longer, self-intersecting loops. 
At /3c, these loops start to percolate through the lattice. For [3 below /3c, the system is described by 
loops of all sizes up to . For all couplings /3, there is no noticeable interaction between the vortex 
loops. 

The correlation functions Cs and Cj can be extended to other objects which are invariant under 
the symmetries of the system. More specific, it would be possible to extend Cs and Cj to measure the 
behavior of monopoles in three (or more) dimensional U(l) gauge theory. The field theoretic connection 
of Cs would then allow us to study the behavior of monopoles in these theories, hopefully including the 
mechanism of monopole condensation. The more obvious use of Cs and Cj for Dirac strings will fail, 
as these strings are gauge dependent. 
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Figure 1: Spin configuration as given in (5). Tlie spins are in tlie 1-3 plane. The links j following from 
these spins are shown in gray. A link jfi(x) = 1 is chosen to point away from x. 



x<^I>y 

Fig. 2a 




Fig. 2b 



Figure 2: Possible diagrams that constitute < ip'^ (x)ip(x)ip'^ (y)ip{y) >. Upper, figure 2a, a diagram 
leading to & = 2 in equation 7c. Lower, figure 2b, a diagram that describes meson-like behavior of the 
aforementioned correlation function, and would lead to & = 1. 
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Figure 4: The correlation functions Cs(r) and versus r, from a 16^ lattice for /3 = 0.30 (figure 4a, 
upper), and for /3 = 0.55 (figure 4b, lower). 
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Figure 5: The exponent a from different lattice sizes, depicted versus /3. 
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Figure 6: The constant c from different lattice sizes, depicted versus /3. in the range /3 = 0.00 to /3 = 0.60 
(figure 6a, upper), and for /3 in the range /3 = 0.35 to /3 = 0.60 (figure 6b, lower). 
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Figure 7: The results for the loop density SI = i^/N^ on a linear scale (top left, figure 7a) and on a 
logarithmic scale (top right, figure 7b), and the average loop length A (bottom left, figure 7c) from all 
four lattice sizes, and the fractal dimension rfj from 16^ only (bottom right, figure 7d). All quantities 
are shown versus /3. 
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Figure 8: The loop density p(l) versus / from 16^ lattices for /3 = 0.30 (upper, figure 8a), and for 
/3 = 0.60 (lower, figure 8b). 
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Figure 9: The behavior of the exponent a (left scale), and the power h (right scale) from fitting p(l) to 
(7c), from 16^. Both parameters are plotted versus /3. 
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Figure 10: The average length of the largest loop, £, versus /3 on a log-plot, from all four lattice sizes. 
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05 
05 


0.48 


0.59l°:f5 


1 


q + o 


7 
5 





008 + 0.005 
UUS_o 005 





221° 


06 
06 


0.49 


o.84l°:i° 


1 


1+0 

-■--0 


3 
3 





QQC+0.003 
'-"-''J-0.003 





251° 


04 
02 


0.50 


O.Qliro? 


1 




3 
2 


0022+° 

U.UUZZ_o 0013 





231° 


03 
03 


0.525 


1 40+0.12 





8t°o 


4 
3 


0005+0 ooO'i 

U.UUUO_o 0003 





28t°o 


05 
05 


0.55 


1 Yq + 0.08 
' "^-0.06 







3 
2 


0001 5+° °°°i5 

U.UUUiO_o 00016 





291° 


03 
03 


0.575 


1 81+006 


0.39lo 


11 
11 






0.288lo 


019 
018 


0.60 




0.01° 


3 
3 






0.351° 


04 
05 


0.65 


q 1O + 0.10 

"^•-'-''-0.08 


0.01° 


3 
3 






0.34t° 


04 
05 


0.70 


4 12+0 16 












0.501° 


09 
09 


0.75 


o 02 + 16 












0.23lo 


04 
05 



Table 1: The results from fitting Cs to (7c) for 6^. 
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a 




h 




c 


d 




0.00 




1 


Q + l 1 





i^n A -1-0 015 

364Io.oi5 


r\ 1 ^-1-0 
O.I6I0 


05 
05 


0.05 




1 


7+0.5 
' -0.4 





C A A -1-0 01 4 

'^^^-0.014 


^ r-T-l-O 

O.lTio 


03 
03 


0.10 




1 


o + O 9 
^-0".4 





ot-*/i-l-0 01 2 
'^^^-0.012 


0.16+0 


03 
03 


0.15 




1 


r\-\-0 9 

^-0.5 





^ ^ r-l-0 Oil 

"^-"-"J-O.Oll 


r\ 1 ^ + 

0.16_o 


04 
04 


0.20 




1 


o + O 5 
^-0.3 





„_„ + 007 
^ ' ^-0.008 


r\ 1 T+O 

0.17Io 


02 
02 


0.25 




1 


r\-l-0 7 

^-0.4 





i-> - n-l-O 01 1 


0.17lo 


03 
03 


0.30 




1 


O + O 4 

f'-o.s 





^ r\-l-0 007 
^-'-'-'-0.007 


0.17+0 


02 
02 


0.35 


r\ r\ A -t-0 12 


1 


Q + 5 

f'-o.s 





IRO+O 008 
-"-^^-0.009 


r\ 1 + 

0.19+0 


03 
03 


0.40 


r\ 4 -1-0 1 


1 


r-l-0 5 

"J-o.s 





1 1 o + O. 008 
-'--'-°-0.008 


0.23^0 


02 
03 


0.41 


^ r\-|-n 1 2 

0-10loi5 


1 


7+0.5 
' -0.4 





1 (T^ + O 009 
-'-'-''J-0.010 


0.20to 


03 
03 


0.42 


r\ -1-0 09 

0-21Io:o7 


1 


c+0.3 

"J-o.s 





oor, + 0.007 

^0^-0.007 


0.23^0 


03 
03 


0.43 


0.221^;^= 


1 


4+0.3 
^-0.2 





nRO + O 005 
'-'^^-0.005 


0.23^0 


02 
03 


0.44 


t-if-,-1-0 1 6 

0-28l°ii 


1 


q + 0. 6 
"^-0.3 





O.q + 0.008 

'-'^"^-0.009 


^ ^ , _i_n 

0.24Io 


05 
06 


0.45 


0.32l°23 


1 


q + 0. 9 
"^-0.5 





025+° °" 
'-'^'J-o.oio 


0.24^0 


07 
08 


0.46 


or'-l-O 18 


1 


q + 0. 7 

"^-0.4 





014+0 007 
'-'-^^-0.007 


0.23^0 


05 
07 


0.47 


^ rr\-l-0 12 


1 


9 + 0.4 
^-0.2 





OOt^ + 003 
'-'U'J_o.003 


0.23^0 


05 
04 


n 4s 

W .T;0 


u.oa_o 12 


1 


9 + 0.5 
^-0.4 





OOq + 0.003 

'-"-'"^-0.003 


9'^+° 
w . _ 


06 
05 


0.49 


7fi+o 09 

U.<D_o 07 


1 


1+0.3 
-■--0.2 


o.ooo9l°:ro^ 


0.23lo 


04 
03 


0.50 


1 08+° °9 





O + 0.4 
°-0.3 


o.ooo7t°:ro^ 


0.28t^ 


03 
04 


0.525 


1 Q4+0.05 


0.85iri^ 






0.254t° 


018 
018 


0.55 


1 co+O 09 
-'-•'^"^-0.06 


0.9l°o1 






0.23lo 


03 
02 


0.575 


2 29+° °® 


0.2l°i 






0.351° 


04 
04 


0.60 












0.35t°o 


05 
05 


0.65 


3 04+0 07 
"^•'-'^-0.07 










0.30lo 


03 
03 


0.70 












0.200lo 


017 
017 


0.75 


3.80l°o:i° 










0.191° 


02 
02 


0.80 












0.128t^ 


014 
013 



Table 2: The results from fitting Cs to (7c) for 8^. 
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a 




h 


C 


d 




0.00 


— 


1 


(j + 0.4 


n S49+0.007 

U.01:»_o 006 


0.181° 


02 

03 


0.05 




1 


7+0.7 
' -0.4 


0.338l°°{2 


0.17t° 


04 
04 


0.10 




1 


7-t-u. 1 
' -0.4 


n ^1 1 +0.010 

U.;>ii-o.oll 


0.18lo 


04 
03 


0.15 





1 


7+0.7 
' -0.4 


9Qn+o.oio 
u.z»u_o 010 


0.18l° 


03 
04 


0.20 





1 


7+0.7 
' -0.3 


oRq+o.oii 
u.zuo_o 010 


0.18lo 


04 
04 


0.25 




1 


7+0.4 
' -0.3 


n oqi +0.007 

U.zoi_o 007 


0.18t^ 


03 
02 


0.30 





1 


7+0.5 
' -0.2 


n 1 qq+0.006 

U.l»»_o 007 


0.18lo 


03 

02 


0.35 




1 


7+0.5 
' -0.3 


n 1 58 + 0.008 
U.idO_o 007 


0.18l° 


03 
02 


0.40 


^ . - 1 n n g 

n 1 o+u uo 
u-J-^-o.os 


1 


(j + 0.2 


^ . ... 1 n n n Q 
n 1 1 t^ + 0.003 
U.iid_o 005 


0.2lt° 


02 

02 


0.41 


n 1 q + 0.09 
'-'•-'-"^-0.08 


1 


r, + 0.4 


0.10lt°oZ 


0.2lt° 


03 

03 


0.42 


n 1 7+0 06 

U-J- ' -0.05 


1 


c + 0.2 
"^-0.2 


0.086l°:°°t 


0.22t° 


02 

02 


0.43 


1 R + 


1 


4+0.3 
^-0.2 


079+0.004 

U.U ( ^_o 004 


0.22t° 


02 
03 


0.44 


o.i6l°:i4 


1 


4+0.8 
^-0.4 


0.056l°;°{{ 


0.22t° 


06 
06 


0.45 


o-2ot°:J^ 


1 


q + 0.4 
"^-0.3 


0.030t°;°°^ 


0.23t55 


04 
04 


0.46 


0.25l°:lo 


1 


q + 0.4 
"^-0.4 


1 n n n 

01 9+0.006 


0.22t° 


06 
05 


0.47 


o-44l°:i^ 


1 


9 + 0.4 
^-0.3 


1 n n n 1 v 

001 q+0.0017 
u.uuia_o 0018 


0.23lo 


05 
05 


U.48 


n t-q+0.06 

'-'•'J'J-O.OS 


1 


9 + 0.2 
^-0.2 


n nni n+o.ooo7 
0.0010_o^ooo7 


t\ oq+O 

U./Olg 


02 
03 


0.49 


0.66t°:°^ 


' -0.17 


00024+° °°°20 
u.uuuz'i.o 00020 


0.220t^ 


019 
019 


0.50 


63+° °^ 


i.48l°:fe 


00001+° °°°°^ 

U.UUUUi_o 00004 


0.191° 


03 
03 


0.525 


1 iq+006 

-'-•-'-^-0.04 


0.961°.{I 




0.23lo 


03 
03 


0.55 


1 ci+o.oe 

-'-•'^-'--0.06 


91+0.17 




0.23lo 


02 
02 


0.575 


9 ^7+0. 06 
^•U ' -0.05 


0.37l°:{^e 




0.29lo 


03 
03 


0.60 


2 54+0 05 

^•"J^-0.04 


o.ool°:r. 




0.34l° 


03 
02 


0.65 


2 96+° 

^•^^-0.05 








0.311° 


02 
02 


0.70 


q 40 + 0.08 
"^•^"^-0.07 








0.26lo 


03 
03 


0.75 


3.68l°:°t 








0.2041° 


010 
009 


0.80 


4.161°:°? 








0.1801° 


015 
014 



Table 3: The results from fitting Cs to (7c) for 12^. 
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a 


b 






c 




d 




0.00 




1 c + l 


2 
4 





O/in + 0.017 

340I0.017 





18l°o 


07 
06 


0.05 




1 n + O 

l-9lo 


9 
4 





(^^^^_|_0 010 
"^^^-0.011 





18t°o 


05 
05 


0.10 




1 r' + O 

1-6-0 


9 
4 





„„r. + oil 

"^'-'"^-0.012 





18l°o 


05 
05 


0.15 




1 ^+0 
1-6-0 


5 
3 





r,Qtr+0.007 
^°'^-0.007 





18l°o 


04 
03 


0.20 




1 /^ + 

1-6-0 


5 
3 





nr'n + 007 
^^^-0.006 





18l°o 


03 
03 


0.25 




1 T + O 

1-7-0 


3 
3 





r,q/l+0.005 

^"5^-0.005 





18l°o 


03 
02 


0.30 




1 T + O 

1-7-0 


6 
4 





iqt; + 0.008 
-^^"^-0.008 





19l°o 


03 
04 


0.35 




1 T + O 

1-7-0 


2 
2 





1 (-4+0.005 
-'-'J^-0.004 





19l°o 


02 
02 


0.40 


0.09+0^04 


1 r/^ + O 

1.56Io 


16 
14 





1 1 r, + 0.002 

-'--'-^-0.002 


0.204^0 


015 
015 


0.41 


0-10-0.05 


1 r+0 

1-5-0 


2 
2 





nqs+0.004 

'-'^-0.003 


0.208lo 


015 
016 


0.42 




1 c+O 

1-5-0 


4 
3 





089+0-004 

'-'°^-0.005 





2ll°o 


03 
03 


0.43 


r\ 1 n + 0.07 

0-12Io.o6 


1 r+O 

1-5-0 


2 
3 





n7n + 0.005 
U<U_o 005 





2ll°o 


03 
02 


0.44 


„ ^ „ + 07 
0-10Io.O7 


1 A+O 

1-4-0 


3 
2 





nrq + 0.005 

'-''J "^-0.005 





2ll°o 


03 
03 


0.45 


0-07Io.io 


1 /l+O 


6 
4 





n9r+0.010 
'-'^'J-0.010 





2ll°o 


06 
06 


0.46 


r\ r|i +0 09 

0-21Io.o7 


1 0+0 
l-3lo 


3 
3 





r)r)tr+0.003 
'-'U'J_o.003 





22to 


04 
04 


0.47 


0-34Io.o5 


1 rin + O 
I.22I0 


19 
15 


o-ooo8l°:S 





2210 


03 
03 


0.48 


'-'•^'J-0.04 


1 3+0 


2 
2 


0-0001210-00012 





2ll°o 


04 
04 


0.49 


0.591°:°^ 


1.251° 


14 
12 


0-000011°:°^^ 





2ll°o 


03 
02 


0.50 


70+0 


1.35to 


15 
12 









2lt°o 


02 
03 


0.525 


1 18+0 04 


0.98lo 


11 
10 






0.233lo 


019 
019 


0.55 


1 70 + 0.01 
' ^-0.02 


0.40lo 


06 
05 






0.30lo 


02 
01 


0.575 


2 06+0 0^ 


0.39lo 


11 
10 






0.288lo 


019 
018 


0.60 


9 04+0.05 
^•"^^-0.05 


0.33t^ 


17 
15 






0.28t^ 


02 
03 


0.65 


2 97+0 04 

' -0.04 










0.297lo 


016 
016 


0.70 


41+0.05 

"^•^-"--0.04 










0.254^0 


015 
016 


0.75 


70+0.08 

' "^-0.06 










0.20lo 


02 
02 


0.80 


4 02+0 05 

^•'-'^-0.05 










0.155t^ 


Oil 
Oil 



Table 4: The results from fitting Cs (left columns) and (right columns) to (7c) for 16^. 
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# = 6 


N = 8 




# 


A 




# 


A 




0.00 


4.48(6) 


58.2(1.1) 


1.2207(14) 


8.37(8) 


68.4(9) 


1.2208(9) 


0.05 


4.48(5) 


56.3(1.1) 


1.2079(12) 


8.57(11) 


65.2(9) 


1.2057(8) 


0.10 


4.66(6) 


51.2(9) 


1.1929(11) 


8.73(8) 


60.6(7) 


1.1912(8) 


0.15 


4.77(7) 


49.1(1.2) 


1.1767(12) 


9.12(8) 


55.6(7) 


1.1760(10) 


0.20 


4.88(6) 


45.0(9) 


1.1615(12) 


9.30(9) 


52.7(7) 


1.1614(8) 


0.25 


5.06(6) 


41.1(7) 


1.1490(12) 


9.76(9) 


47.3(7) 


1.1471(8) 


0.30 


5.18(7) 


38.1(7) 


1.1314(10) 


10.00(10) 


43.3(6) 


1.1333(6) 


0.35 


5.48(8) 


33.0(8) 


1.1167(12) 


10.84(11) 


36.6(5) 


1.1178(7) 


0.40 


6.07(9) 


25.6(6) 


1.0984(12) 


11.91(13) 


29.4(4) 


1.1011(7) 


0.41 


6.16(11) 


24.7(1.0) 


1.0945(17) 


12.20(17) 


27.9(5) 


1.0976(7) 


0.42 


6.63(12) 


23.3(1.3) 


1.0866(12) 


12.91(17) 


25.3(5) 


1.0940(7) 


0.43 


6.76(10) 


19.6(7) 


1.0811(12) 


13.43(18) 


23.1(5) 


1.0885(8) 


0.44 


6.87(9) 


17.5(6) 


1.0774(13) 


14.51(2) 


19.4(5) 


1.0821(10) 


0.45 


7.11(7) 


15.2(4) 


1.0708(13) 


15.4(2) 


16.2(5) 


1.0750(10) 


0.46 


7.06(8) 


13.5(5) 


1.0639(14) 


16.0(2) 


13.6(4) 


1.0677(12) 


0.47 


6.92(13) 


12.6(5) 


1.0609(11) 


16.18(13) 


11.4(2) 


1.0611(8) 


0.48 


6.92(7) 


10.7(4) 


1.0545(17) 


15.87(12) 


10.4(3) 


1.0564(12) 


0.49 


6.80(7) 


9.3(2) 


1.0486(15) 


15.49(10) 


9.01(14) 


1.0492(9) 


0.50 


6.25(8) 


8.45(15) 


1.0432(14) 


14.65(12) 


8.34(13) 


1.0449(9) 


0.525 


5.35(8) 


7.19(10) 


1.0292(11) 


12.35(14) 


6.97(7) 


1.0343(8) 


0.55 


4.37(8) 


6.41(8) 


1.0229(13) 


10.14(11) 


6.37(6) 


1.0261(8) 


0.575 


3.26(8) 


5.72(7) 


1.0164(10) 


7.91(14) 


5.91(4) 


1.0196(8) 


0.60 


2.60(6) 


5.21(7) 


1.0105(8) 


6.30(8) 


5.64(5) 


1.0151(7) 


0.65 


1.51(4) 


3.86(9) 


1.0035(6) 


3.69(6) 


4.95(5) 


1.0076(6) 


0.70 


0.89(4) 


2.70(9) 


1.0020(6) 


2.08(4) 


4.15(4) 


1.0051(10) 


0.75 


0.52(3) 


1.81(9) 


1.0011(4) 


1.17(4) 


3.16(7) 


1.0017(4) 


0.80 


0.274(18) 


1.05(7) 


1.0005(3) 


0.66(3) 


2.17(8) 


1.0009(3) 


0.85 


0.165(12) 


0.65(4) 


1.0002(3) 


0.37(2) 


1.33(7) 


1.0003(2) 


0.90 


0.119(13) 


0.45(4) 


l.O(-) 


0.232(18) 


0.85(7) 


1.00014(15) 


0.95 


0.053(7) 


0.21(3) 


l.O(-) 


0.120(16) 


0.42(5) 


l.O(-) 


1.00 


0.035(7) 


0.14(3) 


l.O(-) 


0.068(12) 


0.23(4) 


l.O(-) 




Table 5: 


The results 


for #, A, and rfj for N = 


6 and N = 


: 8. 



21 









# = 


12 








N = 


: 16 






/? 


# 


A 




# 


A 


di 


0.00 


23 89 


'1 6"! 


75 2 


'6^1 


1.2208 


'4) 


54.2(2) 


76.9 


:4) 


1.2209 


'3^1 


0.05 


24.61 


'1 6"! 


70 6 


'5^ 


1 9065 


'5^ 


55.7(2) 


72.4 


:4) 


1 9065 


'3^ 


0.10 


25 97 




64 8 


'5^ 


11919 


'5^ 


58.6(3) 


66.3 


:3) 


1.1912 


'3^ 


0.15 








) 


1.1 t uo 




60.4(2) 


61.8 


:3) 


1 1 771 


/ 


0.20 


97 7^ 




55-5 


/ 






63.2(2) 


56.4 


:3) 


1 1 ^9^1 


/ 


0.25 


9Q 97 


'1 6"! 


4Q 7 


'4) 


1 1479 


'4) 


66.4(3) 


50.9 


:2) 


1.1484 


) 


0.30 






44 

4:4: . U 




1 . lOOO 




70.0(3) 


45.3 


:2) 


1 1 '^'^'^ 


) 


0.35 


32.8 


:3) 


37.9 


'3) 


1.1176 


'4) 


74.8(3) 


39.02 


;i8) 


1.1184 


'3) 


0.40 


35.8 


:2) 


31.2 


'3) 


1.1022 


'4) 


82.0(5) 


31.82 


;i6) 


1.1027 


'3) 


0.41 


36.6 


:3) 


29.5 


'3) 


1.0984 


'4) 


84.3(4) 


29.91 


;i5) 


1.0990 


'3) 


0.42 


37.8 


:3) 


97 




1 nQ49 




86.2(3) 


28.38 


15) 


1 nQfi9 




0.43 


39.6 


:4) 


25 3 


'3^ 


1 091 6 

1 .\j tj lyj 


'4) 


89.6(5) 


26.1 


:2) 


1 091 8 

1 .\J tJ lO 


'3^ 


0.44 


41.2 


:6) 


23.2 


'5^ 


1 0865 


'4) 


94.0(6) 


23.7 


:2) 


1 0879 


'3^ 


0.45 


45.2 


:5) 


19.3 


^3) 


1.0813 


^5) 


103.8(1.0) 


19.7 


:3) 


1.0821 


1) 


0.46 


49.5 


:6) 


15.5 


;4) 


1.0742 


;6) 


117.4(9) 


15.1 


:2) 


1.0740 


1) 


0.47 


53.3 


:3) 


11.7 


'3) 


1.0638 


7) 


124.3(5) 


12.32 


;i2) 


1.0670 


1) 


0.48 


53.0 


:3) 


10.55 


:i4) 


1.0599 


;6) 


125.3(4) 


10.61(8) 


1.0603 


1) 


0.49 


51.7 


:2) 


9.24 


;8) 


1.0526 


;6) 


122.6(3) 


9.48 


:6) 


1.0544 


;3) 


0.50 


49.6 


:2) 


8.52 


;2) 


1.0488 


;5) 


117.4(4) 


8.53 


:4) 


1.0488 


;3) 


0.525 


42.2 


:3) 
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Table 6: The results for A, and for # = 12 and N = 16. 
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Table 7: The results for fitting p(l) to (7c). 
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Table 8: The results for fitting p(l) to (7c). 
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looU^bUj 


n /I /I 
U.44 




OVO/ 1 o\ 


Q 1 1 / 1 0\ 


1 vnn/ r;n\ 
1 /UU^oUj 


n /I 
U.4o 


1 OQ / 1 




VI n/ r;n\ 
/ lU^oUj 


1 /I (^n/ 1 r;n\ 
14bU^loUj 


n /I 
U.4d 


1 OQ / 1 /I \ 

izo^i4 j 


o 1 n/ '?n\ 


r: Qfi/ nn\ 

oou^yuj 


1 n/1 n/ onn\ 
lU4U^zUUj 




1 1 o/i 7A 
1 lol^l ( J 


1 7C\( AC\\ 
1 ( Ul^4U J 


ooui^yu j 


1^1 n/1 vnA 


0.48 


95(16) 


140(30) 


340(80) 


390(50) 


0.49 


80(16) 


120(30) 


230(60) 


230(30) 


0.50 


65(11) 


90(20) 


170(40) 


180(40) 


0.525 


45(9) 


62(19) 


90(20) 


91(11) 


0.55 


36(5) 


40(11) 


51(6) 


61(9) 


0.575 


26(5) 


30(10) 


46(10) 


48(6) 


0.60 


21(3) 


27(6) 


33(12) 


43(11) 


0.65 


14(3) 


20(5) 


23(4) 


22.4(1.5) 


0.70 


11(3) 


13.2(1.6) 


16(3) 


17.6(1.2) 


0.75 


10(2) 


11.8(2.0) 


13(2) 


17(3) 


0.80 


7.6(1.2) 


8.8(1.3) 


10.4(1.2) 


13(2) 


0.85 


7.2(1.8) 


7.2(1.6) 


9.4(1.3) 


12(3) 


0.90 


6.2(1.4) 


6.8(1.0) 


8.6(0.9) 


9.6(1.7) 


0.95 


5.2(1.6) 


6.4(1.5) 


7.4(0.9) 


9.0(1.8) 


1.00 


5.0(1.6) 


5.2(1.3) 


7.0(1.0) 


8.4(1.2) 



Table 9: The results for the average of the largest loop, £. 
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